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Abstract 

The unsteady aerodynamic effects due to the separated flow around slender delta wings 
in motion have been analyzed using an extension of the Brown and Michael model, as 
first proposed by Arena. By combining the unsteady flow field solution with the rigid 
body Euler equations of motion, self-induced wing rock motion is simulated. The 
aerodynamic model successfully captures the qualitative characteristics of wing rock 
observed in experiments. For the one degree of freedom in roll case, the model is used to 
look into the mechanisms of wing rock and to investigate the effects of various 
parameters, like angle of attack, yaw angle, displacement of the separation point and 
wing inertia. To investigate the roll and yaw coupling for the delta wing, an additional 
degree of freedom is added. However, no limit cycle was observed in the two degree of 
freedom case. Nonetheless, the model can be used to apply various control laws to 
actively control wing rock using, for example, the displacement of the leading edge 
vortex separation point by inboard spanwise blowing. 
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introduction 

At high angles of attack the leeward flow field for slender delta wings is dominated by a 
highly organized vortical flow structure emanating from the sharp leading edges. The 
vortex sheet shed from the leading edge rolls up into a pair of strong vortices. As the 
angle of attack is increased, these leading edge vortices interact with each other and the 
wing itself, to create a sustained and large amplitude rigid body oscillation called wing 
rock. Such oscillations lead to a significant loss in lift and can present a serious safety 
problem during maneuvers, such as in landing. The maneuvering envelope of an aircraft 
exhibiting this behavior is also seriously restricted because the maximum incidence angle 
is often limited by the onset of wing rock before the occurrence of stalL 
In recent years, the effectiveness of active control to alleviate this problem has been 
explored at the Stanford low speed wind tunneL Wongl 18 ! has demonstrated that one 
degree of freedom wing rock can be suppressed by using tangential blowing as a roll 
control actuator and Pedreirof 15 ! is examining the roll and yaw coupling for two degree 
of freedom wing rock and its e liminati on However, the development of efficient control 
techniques and algorithms is constrained by the limited understanding of the basic 
aerodynamic mechanism and effect of various parameters on this phenomenon. This 
information must be understood so that wing rock may be avoided through design or 
efficiently controlled by active means. Also, due to the various conditions involved in the 
wing rock experiments, one lacks a general aerodynamic model where control laws can 
be tested for various configurations before actual implementation. These were the 
motivations for developing the following aerodynamic model. 


Aerodynamic Model 
Assumptions 

An aerodynamic model is needed to obtain a fast estimate of the velocity and pressure 
fields around a delta wing which when combined with the equations of motion will 
capture wing rock. The vortex model for the vortical flow around the delta wing, first 
suggested by Brown and Michael^ and later expanded to the unsteady case by Arena^l 
is adequate for this purpose and is chosen to investigate the characteristics of wing rock. 
The model developed is an extension of the Arena model in the sense that for the one 
degree of freedom case parameters such as sideslip and displacement of the separation 
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)int were added and yaw motion in the plane of the wing i.e., an additional degree of 
eedom was included. 

fhe actual flow field around the wing is illustrated in figure 1(a). In the model, separated 
flow on the delta wing is represented by a pair of line vortices connected to the leading 
edge by a straight vortex feeding sheet, as shown in figure 1(b). It has been shown 
experimentally! 17 ), that most of the axial vorticity of the leading edge vortex is confined 
to a viscous core region having a diameter of the order of 5 % of the local semi-span. This 
fact justifies using a model where all the leading edge vortex vorticity is concentrated 
into two single line vortices. 

Usually aircraft operate in the range of high Reynolds number, where the viscous effects 
are confined to very thin boundary layers along the surfaces and free shear layers in the 
fluid. Thus, it will be assumed that the only role of viscosity is to provide the mechanism 
for flow separation. It will also be understood that the Mach numbers to be used are 
sufficiently low to assume incompressible flow. 

The present model does not predict vortex breakdown and no attempt is made to include 
this phenomena. Wing rock is observed at angles of attack where vortex breakdown does 
not occur! 2 ) and therefore the dynamic simulations will be applied at incidence angles 
where vortex breakdown does not affect the aerodynamics of the wing. The steep 
pressure gradient between the minimum pressure and the primary separation line causes 
an additional flow separation, which usually takes the form of a small secondary vortex. 
The effect of the secondary vortex on wing rock is small and will not be considered. It 
will also be assumed that the flow field is conical. The conical assumption requires that 
the wing geometry be conical and therefore all physical quantities are constant along rays 
emanating from the wing vertex. For a finite delta wing, subsonic conicality is an 
approximation that stems from ignoring the singular nature of the apex and the trailing 
edge effects. Nonetheless, it has been observed that the subsonic flow past a delta wing is 
approximately conical in a region downstream of the apex and upstream of the trailing 
edge. Later a slender body assumption will be used to simplify the governing equations. 

Wing Geometry and Coordinate System 

The unsteady aerodynamics of an aircraft maneuvering at high angles of attack are very 
configuration dependent. Therefore, to provide some insight into the unsteady 
aerodynamics, we will look into a simple configuration to effectively eliminate any 
configuration effects, namely that of a plain slender delta wing. 
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The schematics of the delta wing and coordinate system to be used in the dynamic 
simulations are shown in figure 2. It is assumed that the wing has zero thickness and is 
mounted on a pivot 

Two coordinate systems are used: one inertial ground fixed frame of reference and one 
moving frame attached to the wing. All measurement and operations made with respect 
to the inertial frame are denoted with capital letters while lower case letters are used for 
the moving frame. The numerical problem is posed in the body fixed coordinate system. 
Therefore the relationship between variables in these two frames must be examined. 


d>(/?,r) = <D(r,r) p(R,t) = p{r,t) 


v,=v r V’=V? 


d_ 

dt 


d_ 

dt 


-(Oxr)V, 


In the inertial frame of reference, the X axis is aligned with the freestream velocity. 
Before the dynamic simulation the wing is moved to its initial position by the following 
sequence: 

(a) a yaw-like rotation around the original Z axis through the angle (5 followed by 

(b) a pitch like rotation around the new position of Y axis through the angle a followed 
by 

(c) a roll like rotation around the final position of the X axis through the angle y. 

Once the initial position of the wing is fixed the dynamic simulation uses Euler angles £ , 
t| and <J> to provide the wing motion. 

The base vectors in the inertial frame and the freestream velocity are related to the body 
fixed frame as follows: 


i 

j 


7 " 

1 


U 

V 

=[<w, 

V/ 

0 
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K 


w 

b 
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where 



cosacos/J 

sin y sin a cos fi- cos y sin fi 
cos y sin a cos /J +sin y sin 0 


cosasin/3 -sin a 

sin y sin a sin /J+ cosy cos P sin y cos a 
cos y sin a sin p -sin y cos P cosy cos a 
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cosrjcos^ 

sin 0 sin 77 cos £ - cos <f> sin £ 
cos #sin 77 cos£ +sin <p sin £ 


cosTjsinf 

sin <p sin 77 sin £ + cos <p cos £ 
cos <p sin 77 sin £ -sin <f> cos £ 


-sin 77 
sin 0 cos 77 
cos <j> cos 77 


Complex Potential 

The Prandtl-Glauert equation 


( 1 -A/l) 


<? 2 <D 


<2r 




<? 2 <D 


2 + ^- 2 + A 2 


= 0 


is valid for supersonic and well as subsonic Mach numbers. If the wing is slender and the 
crossflow effect is dominant, the x derivative must be smaller than the other terms and 
the Laplace equation for the crossflow is obtained. 


, d 2 <t> n 

dz 2 ~ 

Unlike the original Brown and Michael model, due to the asymmetry of the flow field, 
the delta wing cross plane will be transformed by conformal mapping to a circle plane 
and the circle theorem, which allows one immediately to write the complex potential in 
terms of the vortex system and its image will be used. The advantage of this approach is 
that the boundary conditions on the wing surface are satisfied exactly and the time 
dependency can be introduced through the boundary conditions. 

The conformal mapping function is given as 


£=-( <j+V<r 2 -4a 2 ) 

where \ represents the circle plane while a represents the physical plane. Figure 3 . is a 
sketch of the approximated flowfield in the crossflow physical plane and the circle plane. 
For a conical flow, the vortex strength increases linearly in the chordwise direction and 
therefore a feeding sheet is necessary in the model. This sheet is represented by a 
mathematical branch cut so that the potential function is single valued and is represented 
by the dotted line in figure 3 . The cross flow velocities v b and w b are functions of the 
angle of attack, sideslip and roll angles and are defined previously. The steady flow field 
in the circle plane is represented by the superposition of a doublet with flow coming from 
two directions and two vortices. To satisfy the tangency condition, image vortices of the 
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opposite strength placed at a 2 /^ must be used. Using these elements, the steady 
complex potential can be written as 


= 4>+iY = Ww/4~1+v1{+41+'^I»(«-4i)-<T-I“ 

V W V $ ) 



-i^ln(|-^)+/^ln 

ZlXr ZK 



The expression for the velocity at any point can be obtained by differentiating the 
previous equation with !;. 


KuaJyiZ) = 17 = V_/H, = -w/l + TT |+ V i l ~JT 

s J V S 

. r 1 . r r l 

-i— — — +i — , 

2nt;-t; r 2 k c a 

S T 

In order to allow for unsteady motion of the wing the complex potential and the 
boundary conditions must be modified. The governing equation is the same for the 
unsteady case. However, since the Laplace equation does not explicitly depend on time, 
the boundary conditions must be time dependent and should be solved at each time step. 
The tangency condition for the unsteady case should be modified to state that the local 
fluid velocity normal to the wing should be equal to the local velocity of the wing itself. 
To satisfy this condition in the unsteady case, a potential function must be superimposed 
with the steady case to account for the unsteady boundary conditions. The derivation for 
this unsteady velocity potential is similar to that of Bisplinghoff, Ashley and Halfman^l 
for the unsteady flow on a pitching airfoil. The unsteady condition due to roll of the wing 
can be satisfied by using a source-sink sheet on the circle and the unsteady condition due 
to yaw can be satisfied using a doublet and freestream. The complex potential for a 
source-sink sheet around the circle and the doublet with freestream from the yaw can be 
expressed as 

Fun,*** = ^f 0 *H+ln(i;-ae ie )add-(x-M:r)n t 

The unsteady flow field can be obtained by differentiating the expression with respect to 

& 
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w„ 


unsteady ' 


= — f* -(*- Z-Cr)n( 1 - 1 

2/r J <> £-««'* *1 
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2 

2 
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The unknown source-sink strength per unit length H + is given in terms of the local 
velocity of the wing as H + =2V W . The expressions will be transformed to the circle plane 
where 

V -U^ = ^ =2aQ x cos0 

= 2yQ I sin0 = 4oQ x sin 0 cos 0 

The expression for the unsteady complex velocity in the circle plane can be written as 

A- 2 / 


w„ 


unsteady 
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C^^e-( x -x Cr) ci, 

f a 2 ) 

l-JT 

Jo Q - ae 

l $ J 


Therefore, the total complex velocity in the circle plane can be written as 


f a 2 ^ 


( a 2 \ 


+ V fc 
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{ ?) 


.r, 1 .r 1 

+i—~ — z — i — 


— i- 
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+i n 


1 4 a 2 Q. f 2 * sin 0 cos 0 


2n . tr 
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Jo 


£ - ae" 


-f-de-{x-ACr)n : 


s' 
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To obtain the pressure distribution on the wing, the tangential velocity due to the 
unsteady motion must be known. However, the equation for the unsteady velocities is 
singular on the wing surface and cannot be used. Therefore, the expression for the 
tangential velocity on the circle plane due to the roll and yaw given by Bisplinghoff et 
al.[3] will be used 


V = — [ 

' irJC 


2r* V. sin 2 0 


-di} + 2(x - A Cr)Cl t sin 6 


n Jo cosd-cosd 

By substituting the expression for the wing velocity in the physical plane, the tangential 
velocity can be found. 


V = 


4Q,a r* cos # sin 2 1 ? 


'i 


t/tf + 2(x - A Cr)Cl z sin 6 = — ^ J c 


Q r a r* cos# -cos # 


n Jo cos#-cos0 ' z n Jo cos#-co sd 

+2(x - A Cr)Q t sin 6 = -2Q,a(l - 2sin 2 0 ) + 2(x - A Cr)Q t sin 6 
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Kutta Condition 

Normally, the separation point for a slender sharp edged delta wing is fixed at the leading 
edge. However, to study the effects of displaced separation, it is assumed that the 
separation point is located at an arbitrary point <r s . This requirement of a separation point 
at <r s is expressed in the present formulation as 


W(a,)=w(i)^- 

do a 



= 0 


Therefore, for the left and right separation points we obtain the following two conditions. 


-iw b 


f, 

1 + T 


f a 2 ) 

+ v b 

1-— y 

£ 2 J 

0 



, .r, 1 .r, 1 .r, i 
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+! 


: r r 1 


r u sin 0 cos 0 


2 n - a 

g, ~z 


n 


-i 

k 




*-fde-{x-u:r)a z 


' a 1 ^ 
1 — ^ 

S. 2 J 


= 0 


Note that the integral is singular but has a finite value at = ae ‘ 9 . Hence to avoid any 
numerical problems an adaptive integration scheme is used. 


Force Free Condition 

In an actual fluid, the fluid pressure is continuous everywhere. However, for the present 
model the feeding sheet is modeled by a branch cut which gives a jump in the potential 
function which in turn creates a pressure discontinuity. The force-free condition is based 
on the fact that there can be no unbalanced forces in the fluid, hence the pressure forces 
from the feeding sheet must be balanced by the force on the vortices. 

The pressure jump across the feeding sheet in terms of the pressure coefficient can be 
shown to be 



The force exerted by the feeding sheet can be obtained by integrating AC p from the 
separation point to the vortex and multiplying by the dynamic pressure. 


r Uft feeding sheet 




The force exerted on the vortex itself is given by the Kutta- Joukowski theorem 
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left vortex 


= ipT, 


W{c,)-u t 


do, 

do, 

dx 

dt 


^ right vortex "" *P^r 




*J 


where 


do. 


dt 


do. 


is the velocity of the vortex with respect to the inertial frame and u b ~r L is 

dx 


due to the inclination of the line vortices with respect to the freestream velocity. 

From the force free conditions, the force from the feeding sheet must cancel the force 
from the vortex. Combining the equations for both the left and right vortex we have 


do. 


dt 


= -u. 


do, 

dx 


•+* '(ff.)-k 


dx dt 


I o-o. 


To simplify the equations, the conical flow assumption will be used. For conical flow the 
vortex strength and position are linearly increasing functions of x. Therefore, 


do _ o,\ 

— - = — tane 


dr v r v . 

— - = — tane 
dx 2 a 


dx 2 a 

All the equations up to this point have been defined with respect to the ground fixed 
inertial frame. Since the wing will go into motion, it is convenient to transform the 
equations to the body fixed frame. We have previously shown that for the two ground 
fixed and body fixed coordinate systems 


do. 


dt 


do. 


dt 


+Qxr 


All equations will be solved in the circle plane and can be transformed by 

.2 




a_ 


do. 

1 

* 

h 

Z~a 2 dL 

dt 

r dt 

r SI ^ 


Milne-Thompsonl 11 ! has shown that one must be careful when transforming equation 
carrying vortex expressions to different planes, because the solution will not be the same 
when transformed back to the physical plane. Therefore to calculate the velocity at the 
vortex center, the following terra must be added. 




e ±i r is 




The final form of the force free condition can be written as 
_ £ £ 


dl 

dt 


Z-a 2 


w(0- 


JL. 


u b tane 

2nl 


la 


2 |&+1 Z.+-T 

W l £ 
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r 17 ( -2\T| EiP+a 1 ) 

— “ £ + — - 4 , + — t — -r--(j:-ACr)i2 

r,[{ L) r (,) JJ (^-a 1 ) 

The final step is to transform the equations into conical variables to eliminate the 
chordwise dependence of the steady terms. The vortex position, strength and time will be 
non-dimensionalized by the following variables. 

£ r rV_ 

a 2naV Cr 

The final system of equations is 


-A 1 + ‘ W l-±U 


v -i a u a c-ff c ' c-e 


2£tanrfll f«sin9cos£ rf9 _|' JL-Aj; 1--L U 0 

*Cr J» Z-f ICr ) 'I £*J 


£ «ste)#T*iAr ^-54 4e*±Hc4 


dt \ r z : 1 - 1 v v (^-i ) 2 xtan£ v - * 


ir.Tft. O f=. 1 'll! £fe +1 ) 2 (, ACr^. 


where 


w ( f -)= -iiif i+J- l+i-l _& s ^— + K '~ 

>_( rj v-l rJ f-ff r-c r .j_ 

k. €r 


2xtan£ft' r2*sin0cos0 , Q ( x -Y-.f, 1 ^ 

J. {« r{ l ~Y) 


We now have six unknowns (two normal and spanwise vortex positions and two vortex 
strengths) and six equations (two from the Kutta condition and four from the zero force 
condition). 


page 17 



Equations of Motion 

For a thin delta wing mounted on a ball and socket sting at point c, as shown in figure 2, 
the Euler rigid body rotational equations of motion are of the form 

I c =/ ( n+Qx(/ c Q) 

where M c is the external moments around point c, £2 is the wing angular velocity vector 
and I c is the inertia matrix about c. For the special case of principal axes. 


where 


/< = 


u 
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0 0 
K o 
o /« 


7 „ = 7 P»Cr* tan 3 e I„= p w Cr 4 tan e 


18 
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7 a =I xx +I yy = 


tan 2 e 




pjCr* tan e 


The equation of motion can be written in component form 

M x = I xx^x + ( 7 « — ^yy)^Py 


M y ~ I yy^y ( 7 « 7 n 
A/, + (/„-/ Jf2,£2, 

The components of £2 can be written in terms of Euler angles using the following 
relations. 


0 / 


-sin 77 

0 
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c" 

n. 

- 

sin 0 cos 77 

COS0 

0 

J7 

Pz. 


cos 0 cos 77 

-sin0 

0 

0 


Differentiating these equations and noting that I zz =I yy -i-I xx . the equations of motion can 
be rewritten as 

M x = I„ (0 - sin tj£) + /„ [cos 0 sin 0(cos 2 - 17 2 ) - 2 cos 77 sin 2 07 )£ j 

M y = I n (cos 077 +sin 0 cos 77 £)+/ w [sin r/(cos 17 cos 0 £ 2 - 2 sin j 

M t = I a (cos 0 cos rj£ - sin 0 *))+ 7 «( sin 77 cos 77 sin 0£ 2 -2 sin 0 cos 770£ - 2 cos # 77 ) 

- I n (sin r) cos 1 7 sin 0£ 2 + 2 cos 0 sin 
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The external moment terms are generated by gravity, frictional damping and the 
aerodynamic forces. 


gravity terms: p w Cr 3 tan e ^ - A j#(cos <p cos ry - sin 0 cos T)k ) 

damping terms: -p^-sin rj£)i - p y (cos 07) + sin 0 cos T]^)j - p t [cos<j>cosTj^ - sin(f>rj)k 
aerodynamic moment terms: ^pV 3 Cr 3 tan e(C x i +C y j + C t k ) 


where 


C = 


tane 


rfe - c J sin6cos ^ - c J sinft/0 


c l = o 


c,=l- 


V 1 


The pressure coefficient is obtained using the unsteady Bernoulli equation. 

2U& A, 

Vl K. 2 Vj 

To transform the Cp expression to the body fixed frame, we have previously shown that 

<*D| 


db 

dt 


dt 


+ (Qxr)« V/D 


All the dynamic simulation results which will be shown later do not contain the frictional 
damping term. Nonetheless, these terms can be added to simulate the friction from the 
bearings which is present in the actual experiments. 

The equations of motion can be rewritten in matrix form 


T 


'1 

sin 0 tan 77 

cos0tan?7 

X' 


= 

0 

COS0 

-sin0 


f. 


0 

sin 0/ cos 77 

cos 0/ cos 77 

X. 


where 


R x = — pVjCr 3 —~—C x - — (0 - sin r/£) - cos 0 sin 0(cos 2 rj£ 2 - V 2 ) + 2 cos 77 sin 2 <pi)£ 
2 /_ I __ 

JOC XX 

R y = - pV*Cr 3 ^^-C y (cos077+sin0cos77£)+— Cr J tan£| — -A Wos^cosrj 

- sin rj(cos 77 cos 0£ 2 - 2 sin 
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/?, =^pV'iCr 3 -^^C { - — (cos 0 cos sin ^77)-— Cr 3 tan£ --A^sin^cosTj 

2 ! a ! zz f a U yT 


+ 2— (sin 0 cos 770£ + cos 007?) + — — —sin 77 cos 7/ sin 0£ 2 + 2— cos 0 sin 7777 £ 

For the two degree of freedom case, £ = £ = £ = 0. After some manipulation of the 
equations and changing into non-dimensional form for roll and yaw we obtain 


v. 1 _ 5 tane _ Cr u . . 

<P = 2 PCr i~ Cl ~~I < * > ~ cos< f >sm ^ 




?• ~HiCr i. P w Cr 5 (2 A I • 

t -y 2L -p-tan^--Ajgtan0+2-^-tan00C 

For the one degree of freedom case there is an additional constraint, £ = £ = £ = 0. In 
non-dimensional form, the equation of motion for roll is given as 





tane 


Cr p. 


¥ 


To begin the dynamic simulation, the wing is first set at an initial position. The static 
flow field is then solved and the wing is released. The pressure distribution is calculated 
and the aerodynamic moments are obtained by integrating the pressure distributions. The 
new wing position and angular velocities are found from solving the equations of motion. 
Because the method can only predict unit time steps, a predictor-corrector method by 
Carnahan! 5 ! is used for time marching instead of a Runge-Kutta method. Using the new 
wing position and velocity, the unsteady vortex position and strength are obtained using 
the Kutta and zero-force conditions. A flow chart of the described process can be seen in 
figure 4. 


p«ge20 


Results and Discussion 

For analyzing wing rock aerodynamics, an extremely useful tool is the concept of energy 
exchange proposed by Nguyen et. all 14 l. During wing rock motion, the energy added or 
extracted from the system for a certain time interval can be expressed as 

A £ = ^ pVjCrj' 3 C , (t)mdt = | pVlCrj* C, ( Qdf . 

Therefore, for a steadyjitate wing rock limit cycle the net energy exchange is given by 



From figure 5 it can be seen &at for a counterclockwise loop, energy is added to the 
system therefore the motion is destabilized while for a clockwise loop, energy is 
extracted from the system and the motion is stabilized. The amount of the net energy 
exchange, AE is proportional to the area contained in the loop. 

One Degree of Freedom Wing Rock 

Figure 6 shows the time history of one degree of freedom wing rock where there is no 
sideslip and separation point displacement The angle of attack is 20° and the semi-apex 
angle of the delta wing is 10°. The wing is released at two different initial angles, 0.5° 
and 85°. The wing rock amplitude and frequency is the same for both cases assuring that 
the motion is independent of the initial conditions and is a true limit cycle. 

The steady state total sectional roll moment coefficient plots are shown in figure 7(a). 
Using the energy exchange concept, it is shown that the clockwise destabilizing lobe 
occurs at small roll angles, while the two counterclockwise stabilizing lobes occur at 
large roll angles. The area of the destabilizing lobe is identical to the area of the two 
damping lobes indicating there is no energy exchange and hence steady state is observed. 
By showing the sectional roll moment coefficients for the top and bottom surfaces 
separately in figure 7(b), we can see that all the instability comes from the top surface 
while the bottom surfaces only contributes to the damping. 

The roll moment coefficients plots have shown when and where the damping and 
destabilizing occur but do not tell us the aerodynamic mechanism for this phenomena. 
The physics underlying the motion can be explained by figure 8, the dynamic vortex 
position and strength. Figures 8(a) & (b) show the static and dynamic, spanwise and 
normal position of the vortex. The spanwise vortex position does not vary much from the 
static position during wing rock and leads or lags depending on the roll angle. However, 
the vortex position normal to the wing surface always lags the static position by a large 
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amount, regardless of the roll angle. This vortex position lag gives the destabilizing 
moments which initiate the wing rock motion. The hysteresis is largest at small roll 
angles and this coincides with the roll angles where the destabilizing lobes appear in 
figure 7(a). If there were no damping the wing rock motion would continue to grow 
without bound and therefore damping is necessary to generate limit cycle motion. The 
damping lobes can be explained from the static and dynamic vortex strength shown in 
figure 8(c). The dynamic vortex strength is constantly leading the static vortex strength. 
This leading is largest at large roll angles and this also coincides with the damping lobe 
positions shown in figure 7(a). Steady state oscillation occurs when the instability from 
the vortex normal position lag balances the damping from the vortex strength lead. 

With the mechanisms of wing rock identified, the model was used to investigate the 
effects of various parameters on the aerodynamics of wing rock. 

Effect of Angle of Attack 

It is well known that wing rock does not occur at low angles of attack and as the angle of 
attack is increased an aircraft is more susceptible to wing rock motion. Therefore, the 
effect of angle of attack on the wing rock motion was studied. 

Figure 9 shows that by keeping the wing semi-apex angle the same as in the previous 
case but by lowering the angle of attack to 10°, the wing rock motion will damp and no 
limit cycle will be observed. Figure 10(a) shows the total sectional roll moment 
coefficient for one cycle. The counterclockwise loop direction shows that the whole cycle 
is damped and figure 10(b) shows that both the top and bottom wing surfaces contribute 
to the damping. No instability lobes are present. However, if the angle of attack is raised 
to 30° with the same wing configuration, figure 11 shows that the wing motion will 
diverge. Again by examining the total sectional roll moment coefficient diagram, figure 
12(a), the loop is in the clockwise direction, indicating an energy addition to the system 
and hence destabilization. Figure 12(b) shows that the bottom surface adds damping to 
the motion. However, the top surface instability lobe is so large that it overwhelms the 
bottom surface and gives oscillatory divergence. In actual flight, oscillatory divergence 
will not occur because as the wing rock motion grows vortex bursting will occur. The 
vortex bursting contributes to a sudden loss of lift of the wing and consequently has a 
stabilizing effect on wing rock. The present model does not predict vortex breakdown 
and therefore oscillatory divergence behavior is exhibited. 

The aerodynamic model shows the effect of angle of attack by qualitatively explaining 
why wing rock does not appear at certain low incidence angles and why wing rock is 
more susceptible at higher angles of attack. 
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Effect of Yaw Angle 

To see the effect of yaw on the one degree of freedom wing rock, it is assumed that the 
wing is placed at a certain fixed yaw angle and left free to roll. Even though the wing is 
symmetric, figures 13(a) and (b) show asymmetric amplitudes of wing rock due to the 5° 
sideslip angle of the wing. Again, it can be seen that even for the case with sideslip the 
wing rock amplitude and frequency is independent of the initial roll angle. 

Figure 14(a) shows the total sectional roll moment coefficient during one steady wing 
rock cycle at 5° yaw angle. It is interesting to note that the shape of the curve is almost 
identical to that of figure 7(a), the wing rock case with no yaw. The only difference is 
that the curve in figure 14(a) is shifted to the right by 14°. Again, figure 14(b) gives an 
identical plot to that of figure 7(b) and the only difference is the shift to the right All 
other trends, like the damping and instability lobe sizes or top and bottom: wing surface 
contributions to the wing rock motion are yaw independent It has been shown that the 
wing rock motion is strongly dependent on the angle of attack, yaw angle and semi-apex 
angle of the wing. Therefore, this shifting to the right can be explained by the concept of 
effective angle of attack, effective sideslip and effective semi-apex angle. 

An effective angle is defined as the angle that the wing "feels". For example, even if the 
angle of attack is fixed, as the wing rolls, the angle of attack that the wing feels, i.e. the 
effective angle of attack will reduce as the roll angle increases and will reach 0° as the 
roll angle goes to 90° when there is no yaw. 

From simple trigonometry, the effective angle expressions are given as 

sin a cos cos(y + ft) + sin /?sin(y + ft) 
cos a cos /J 


a 4 = tan 


-1 


w, 


— = tan 
\ u >) 





sin a cos /3(cos y sin ft - sin y cos ft) + sm /J(cos y cos ft + sm y sin ft) 

= tan" — 

coscr cos p 

=£ ~ Peff .right = £ + Ptff 

Figure 15 shows these effective angles for both 0° and 5° sideslip. It can be seen from the 
figure that the effect of the sideslip is to shift the effective angles which determine the 
vortex position and strength. Hence the wing rock motion shifts to the right as the angle 
of sideslip goes from 0° to 5°. The effective angles are also shifted about 14°, which is 
the same amount that the sectional roll moment diagrams were displaced. 
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If we compare the wing rock amplitude for the zero side slip case (figure 6) and the 5° 
sideslip case (figure 13), again we can see that for the positive roll angles the wing rock 
amplitude is increased by an amount of 14° while for the negative roll angles the 
amplitude is decreased by 14°. The dynamic spanwise and normal vortex positions and 
unsteady vortex strength diagrams are shown in figure 16. These show the same shape 
and magnitude as for the zero sideslip case, but are once again shifted to the right by 14°. 
Therefore, the wing rock mechanism for wing rock with sideslip is the same as the zero 
sideslip case, but the effect of sideslip is to shift the roll angles by an amount which 
depends on the amount of the initial yaw angle of the wing. 

Effect of Separation Point Displacement 

It has been shown by Wongf 18 ) that one degree of freedom wing rock can be controlled 
by active means such as changing the leading edge vortex separation point by blowing. 
Although for thin delta wings the separation point is fixed at the leading edge, the model 
was generalized to be able to change the separation point locations and to observe the 
corresponding effect on the wing rock motion. 

The Brown and Michael model shows that the general effect of moving the separation 
point inboard is to move the vortex inward and slightly downward and also to weaken the 
vortex strength. Figure 17 shows the wing rock time history for two initial conditions 
when the right separation point is displaced inboard by 1% of the semi-span. In can be 
seen that there is very little change in the negative roll angle wing rock amplitude 
compared to the no separation displacement case, figure 6, while the positive roll angle 
wing rock amplitude is reduced. The total sectional roll moment diagram, figure 18(a) 
shows that unlike the case of wing rock with sideslip, there is no roll angle shifting. The 
damping lobe areas are identical to the area of the unstability lobe indicating there is no 
energy exchange. However, the enclosed areas of the two damping lobes are different, 
which explains why the wing rock amplitudes for positive and negative roll angles are 
different Once again it can be seen from figure 18(b) that the bottom surface contributes 
only to the damping while all the instability comes from the top surface. For positive roll 
angles, the top surface shows a very small damping lobe while the damping lobe for the 
negative roll angles is quite large. Figure 19(b), the vortex normal position hysteresis 
which gives the unstability at small angles is quite symmetric with respect to 0°. This 
corresponds to the nearly symmetric clockwise instability lobe in figure 18. The vortex 
strength hysteresis responsible for the damping has a small magnitude for positive roll 
angles and a large magnitude for negative roll angles which is again consistent with the 
results in figure 18. Therefore, the asymmetry of the wing rock amplitude for the 
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displacement separation case results mainly from the different size damping lobes caused 
by the asymmetric vortex strength hysteresis. 

The fact that even a small separation point displacement ( 1 % of the serai-span) can give 
a large difference in wing rock amplitude, suggests that moving the separation point by 
inboard spanwise blowing can be a powerful mechanism for controlling the wing rock 
motion. 

Effect of Wing Inertia 

To see the effect of wing inertia, the roll moment of inertia, I** of a wing of semi-apex 
angle 10° at an angle of attack of 20° was varied from 90 kg-m 2 to 210 kg-m 2 in 30 kg- 

m 2 intervals. 

Figure 20(a) shows that the wing rock amplitude is independent of the wing inertia. The 
wing rock frequency was also examined. Figure 20(b) shows that the inverse square root 
of the inertia was almost linearly proportional to the wing rock frequency. This indicates 
that the wing rock behaves like a second order linear oscillator. The same results for the 
inertia variation were seen in the experiments of Arena! 1 ! and once again shows that the 
model captures the qualitative trends seen in experiments. 

Two Degree of Freedom Wing Rock 

Experiments to understand and actively control the two degrees of freedom wing rock in 
roll and yaw are being conducted in the Stanford low speed wind tunnel by Pedreiro. 
However, before an efficient control algorithm to stabilize the motion can be devised, the 
aerodynamics of this phenomena must be fully understood. Also because of the various 
parameters involved in the experiments, a general aerodynamic model were various 
control laws can be tested before actual implementation in the wind tunnel is necessary. 
Therefore, the previous one degree of freedom model was extended to accommodate yaw 
motion in the wing plane as an additional degree of freedom. 

Figures 21(a) and (b) show how the roll and yaw Euler angles and the effective angle of 
attack, effective sideslip and effective roll angle change with time. Unlike the one degree 
of freedom case, no limit cycle was observed. Various initial conditions and damping 
terms of different magnitudes were added to the equations of motion but still no limit 
cycle could be noticed. Figure 22 shows how the top and bottom wing surface sectional 
roll moment coefficients change as a function of roll and yaw angles. The energy 
exchange concept cannot be used, because no organized damping or instability lobes are 
formed. Figures 23 and 24 show how the unsteady spanwise and normal vortex positions 
and strength change with the Euler roll and yaw angles. Unlike the one degree of 
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freedom case, no analysis of the mechanism of the motion can be made by looking at the 
lead or lag of the vortex position or strength. 


pm* 2* 



Conclusions 

Following the work of Arena, an aerodynamic model using an extension of the Brown 
and Michael model has been developed to capture the dominant features of wing rock 
motion. For the one degree of freedom case, the model explains the mechanisms which 
generate and sustain wing rock. Instability which drives the wing rock motion comes 
from the lag in the vortex position normal to the wing surface which mainly occurs at 
small roll angles while the damping comes from the unsteady vortex strength lead. 
Steady state wing rock occurs when the instability from the vortex position lag, balances 
the damping from the vortex strength lead. 

The effect of increasing the angle of attack is to enlarge the instability lobe and to make 
the wing more susceptible to wing rock. The effect of sideslip on the wing rock motion 
has been shown to give a shift in the roll angle and hence give an asymmetric wing rock 
amplitude. The reason behind this roll shifting has been explained by the effective angle 
concept. Moving the separation point also gives asymmetric wing rock amplitudes by 
changing the sizes of the damping lobes. 

For the two degree of freedom case no limit cycle was observed. However, the two 
degree of freedom model can still be used to test various control algorithms. 

Future work can include adding a pitching motion of the wing to create a three degree of 
freedom wing rock model. 
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Circle Plane 


Figure 3. Sketch of Physical and Circle Plane 



Figure 4. Flow Chart for Dynamic Simulation 
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Figure 6(b). Wing Rock Time History (<t> 0 =85°, a=20°, 0=0°. e=10°, 5=0%) 



Figure 7(a). Total Sectional Roll Moment Coefficient (a=20°, (5=0°, e=10°, 5=0%) 
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Figure 7(b). Top and Bottom Sectional Roll Moment Coefficient (a=20°, (3=0°, e=10°, 5 

= 0 %) 
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Figure 8(a). Span wise Vortex Position During Wing Rock (ct=20° f (3=0°, e=10°, 5=0%) 
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Figure 8(c). Unsteady Vortex Strength During Wing Rock (a=20°, ($=0°, e=10°, 5=0%) 
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Figure 14(b). Top and Bottom Sectional Roll Moment Coefficient (a=20°, P=5°, e=10°, 

5=0%) 
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Figure 15. Effective Angle of Attack, Yaw Angle and Semi- Apex Angle as a Function of 

Roll Angle for (J=0° and 0=5° 
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Figure 16(a). Spanwise Vortex Position During Wing Rock (oc=20°, P=5°, e=10°, 6=0%) 



Figure 16(b). Normal Vortex Position During Wing Rock (a=20°, (5=5°, e=10°, 5=0%) 
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Figure 17(a). Wing Rock Time History 
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Figure 18(a). Total Sectional Roll Moment Coefficient (a=20°, (5=0°, e=10°, 8p0%, 8 

r=l%) 
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Figure 18(b). Top and Bottom Sectional Roll Moment Coefficient (a=20°, (5=0°, e=10°, 

5 1= 0%, 8j=l%) 
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Figure 19(a). Span wise Vortex Position During Wing Rock (ot=20°, (S=0°, e=10°, 5p0%, 

5r=l%) 


page 45 









Figure 19(c). Unsteady Vortex Strength During Wing Rock (a=20°, {J=0°, e=10°, 8 

1=0%, 5r=l%) 
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Figure 20(a). Effect of Inertia Variation on Wing Rock Amplitude (oc=20°, |J=0 o , e=10°, 

5=0%) 
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Figure 20(b). Effect of Inertia Variation on Wing Rock Frequency (ot=20°, P=0°, e=10°, 

5=0%) 
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Figure 21(a). Row and Yaw Euler Angle Time History (a=10°, (5=0°, e=5°, 5=0%) 



Figure 21(b). Effective Angle of Attack, Yaw Angle and Roll Angle Time History (ot=10 

°, P=0°. e=5°, 5=0%) 
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Figure 22(a). Top and Bottom Sectional Roll Coefficient as a Function of Roll Angle (a 

=10°, £=0°. e=5°, 5=0%) 



Figure 22(b). Top and Bottom Sectional^ Roll Coefficient as a Function of Yaw Angle (a 

=10°, 0=0°. e=5 v , 5=0%) 
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Figure 23(a). Spanwise Vortex Position as a Function of Roll Angle (a=10°, (3=0°, e=5°, 

8 = 0 %) 



Figure 23(b). Normal Vortex Position as a Function of Roll Angle (a=10°, (3=0°, e=5°, 8 

= 0 %) 
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Figure 24(a). Spanwise Vortex Position as a Function of Yaw Angle (a=10°, P=0°, e=5°, 

5=0%) 
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Figure 24(b). Normal Vortex Position as a Function of Yaw Angle (<x=10°, (J=0°, e=5°, 5 

= 0 %) 



Figure 24(c). Unsteady Vortex Strength as a Function of Yaw Angle (a=10°, (5=0°, e=5°, 

5=0%) 
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